Accepted for publication in ApJ 

Preprint typeset using L£T^]X style cmulateapj v. 5/2/11 



SPECTROSCOPY OF BROAD LINE BLAZARS FROM 1LAC 

Michael S. Shaw 1 , Roger W. Romani 1 , Garret Cotter 2 , Stephen E. Healey 1 , Peter F. Michelson 1 , Anthony C. 
S. Readhead 3 , Joseph L. Richards 3 , Walter Max-Moerbeck 3 , Oliver G. King 3 , William J. Potter 2 

Accepted for publication in ApJ 



(N 

o 

(N 



ABSTRACT 

We report on optical spectroscopy of 165 Flat Spectrum Radio Quasars (FSRQs) in the Fermi 
1LAC sample, which have helped allow a nearly complete study of this population. Fermi FSRQ 
show significant evidence for non-thermal emission even in the optical; the degree depends on the 
7-ray hardness. They also have smaller virial estimates of hole mass than the optical quasar sample. 
This appears to be largely due to a preferred (axial) view of the 7-ray FSRQ and non-isotropic 
(H/R ~ 0.4) distribution of broad-line velocities. Even after correction for this bias, the Fermi FSRQ 
show higher mean Eddington ratios than the optical population. A comparison of optical spectral 
properties with Owens Valley Radio Observatory radio flare activity shows no strong correlation. 
Subject headings: galaxies: active — Gamma rays: galaxies — quasars: general — surveys 
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1. INTRODUCTION 

The Fermi Gamma-Ray Space Telescope was launched 
on 2008 June 11 . Its primary instr ument is the Large 
Area Telescope (jAtwood et all 120091 LAT). Fermi gen- 
erally operates in sky survey mode, observing the entire 
sky every 3 hours, and providing approximately uniform 
sky coverage on time scales of days to years. 

The Fermi LAT First Source Catalog (jAbdo et all 
l2010al 1FGL) catalogs the 1451 most significant sources 
detected in Fermis first year of operation. Based on the 
1FGL catalog, The First C atalog of AGN Detected by 
the Fermi LAT (jAbdo et alll2010cL 1LAC) is the largest 
radio-7 selected sample of blazars to date, associating 
671 7-ray sources to 709 AGN (some may be unresolved 
composites) in the high-latitude sample. 

Our quest is to optically characterize these sources, 
seeking maximum completeness in spectroscopic identi- 
fications and using the spectra to constrain the prop- 
erties of these AGN. Optically, the Fermi sources are 
evenly split between Flat Spectrum Radio Quasars (FS- 
RQs) and BL Lacerate Objects (BL Lacs). In this paper, 
we focus on the FSRQs. A companion paper (Shaw et 
al, in prep.) addresses the BL Lac objects. 

In <j2j we discuss the observational program and the 
data reduction pipeline. In <J2 we describe the measure- 
ments and derived data products. In 21 we measure the 
continuum emission and non-thermal pollution. In fj5[ 
we estimate the black hole masses and Eddington ratio 
of the Fermi FSRQ. In $6l we discuss the orientation 
and shape of the the broad line regions in this popula- 
tion, and in $71 we relate this data set to on-going radio 
monitoring of these AGN. 

In this paper, we assume an approximate concordance 
cosmology — tt m = 0.3, il\ — 0.7, and H = 70 km s _1 
Mpc -1 . 

2. OBSERVATIONS AND DATA REDUCTION 
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2.1. The FSRQ Sample 

This paper reports on a multi-year observing campaign 
to follow-up the Fermi blazars. A principal aim is to 
achieve high redshi ft completeness for the 1LAC sample 
(jAbdo et al.H2010d ). 

In this paper, we discuss the spectra of FSRQs and 
other LAT blazar associations with strong emission lines. 
A major contribution is new spectroscopy of 165 of 
these blazars. To extend the analysis, we also measured 
archival spectra of 64 SDSS blazars in the sample, for a 
total of 229 spectra. 

This work takes the 1LAC high latitude sample to 
96% type completeness, with 316 FSRQs, 322 BL Lacs, 
33 other AGNs, 4 LINERs, and 4 Galaxies. There 
are 30 remaining associated flat spectrum radio sources 
of unknown type - generally these represent objects 
that are optically extremely faint (R > 23) or show 
faint continuum-dominated spectra, where current spec- 
troscopy does not have sufficient S/N to unambiguously 
confirm a BL Lac- type ID. 

The most important sub-set of this emission-line sam- 
ple are the objects with traditional FSRQ properties - in 
addition to the flat spectrum radio core emission which 
allows the LAT counterpart association, we require emis- 
sion lines with kinematic FWHM > 1000 km s _1 and 
bolometric luminosity > 10 42 erg s" 1 . We find that 
188 FSRQ meet these criteria, including 10 low-latitude 
sources with 1LAC FSRQ associations. In addition, some 
11 BL Lacs show well-detected broad lines. For this pa- 
per we adopt the traditional heuristic BL Lac definition: 
continuum-dominated objects with observed frame line 
equivalent width (EW) of < 5 A and, where meas ured, 
Balmer break strength of < 0.5 (jHealev et al.ll2008f ). We 
classify an object as a 'BL Lac' if it meets these spectro- 
scopic criteria at any epoch. For 6 of the 11 BL Lac our 
spectra includes epochs in a 'low' state where decreased 
continuum reveals broad emission lines with > 5 A EW. 
The other 5 objects satisfy the BL Lac criteria in all 
of our spectra, but nevertheless show highly significant, 
albeit low EW, broad lines. 

The emission line sample contains 29 other objects - 
spectroscopically these are 9 galaxies, 5 LINERs, and 15 
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other AGN. These show only strong narrow lines. While 
the line strengths rule out BL Lac IDs for these objects, 
they are manifestly different from our typical FSRQs. 
They may be misaligned radio galaxies, intrinsically weak 
AGN, narrow-line Seyferts and other less common types 
of 7-ray emitters (jAbdo et al.ll2010bh . 

Note that this sample of digital spectra does not in- 
clude measurements for a number of bright, famous 
blazars, with historical spectroscopic classifications in 
the literature. Indeed, 137 of the 316 Fermi FSRQs 
fall into this category. Measured in R, the sources with 
archival spectra are brighter by 1.4 magnitudes than our 
sample, so their omission may introduce systematic ef- 
fects, discussed briefly in Section [531 



2.2. Observations 

We have used medium and large telescopes in both 
hemispheres in a many-faceted assault on this sam- 
ple. Observations were obtained from the Marcario 
Low Resolution Spectrograph (LRS) on the Hobby- 
Eberley Telescope (HET), the Large Cass Spectrom- 
eter (LCS) on the 2.7m at McDonald Observatory, 
on the ESQ Faint Object Spectrograph and Camera 
(iBuzzoni et alJT98l EF OSC2) and ESO Multi-Mode In- 
strument (pekker et all 1 19861 . EMMI) at the New Tech- 
nology Telescope at La Silla Observatory (NTT), on the 
Double Spectrograph (DBSP) on the 200" Hale Tele- 
scope at Mt. Palomar, on the FOcal Reducer an d low dis- 
persion Spectrograph (|Appenzeller et all 119981 . FORS2) 
on the Very Large Telescope at Paranal Observatory 
(VLT), and on the Low Resolution Imaging Spectrograph 
(LRIS) at the W. M. Keck Observatory (WMKO). Ob- 
servational configurations and objects observed are listed 
in TableHJ 

The observing runs jointly targeted emission line ob- 
jects discussed in this paper and BL Lacs, to be discussed 
in a companion paper (Shaw et al, in prep.). All objects 
discussed in this work have highly significant emission 
lines and confirmed spectroscopic redshifts. 

All spectra are taken at the parallactic angle, except 
for LRIS spectra using the atmospheric dispersion cor- 
rector, where we observed in a north-south configuration. 
In a few cases, we rotated the slit angle to minimize con- 
tamination from a nearby star. At least two exposures 
are taken of every target for cosmic ray cleaning. Typical 
exposure times are 2x600 s. 

With the variety of telescope configurations and vary- 
ing observing conditions, the quality of the spectra are 
not uniform: resolutions vary from 4 to 15 A, exposure 
times from 360 s to 2400 s, and telescope diameters from 
2.7 m to 10 m. During this campaign, we generally used 
the minimum exposure required for a reliable redshift, 
rather than exposure to a uniform S/N. This may intro- 
duce selection effects into the sample, discussed briefly 



2.3. Data Reduction Pipeline 

D ata reduction was perfor med with the IRAF pack- 
age (|TodyHl986t IVa ldcs 1986) using standard techniques. 
Data was overscan subtracted, and bias subtracted. 
Dome flats were taken at the beginning of each night, 
the spectral response was removed, and all data frames 
were flat-fielded. Wavelength calibration employed arc 



lamp spectra and was confirmed with checks of night 
sky lines. For these relatively fain t objects, we employed 
an optimal extraction algorithm ((Valdes 1992) to maxi- 
mize the final signal to noise. For HET spectra, care was 
taken to use sky windows very near the longslit target 
position so as to minimize spectroscopic residuals caused 
by fringing in the red, whose removal is precluded by 
the rapidly varying HET pupil. Spectra were visually 
cleaned of residual cosmic ray contamination affecting 
only individual exposures. 

We performed sp ectrophotom etric calibration u sing 
standard stars from lOkd (|1990D and IBohlinl (|2007f) . In 
most cases standard exposures were available from the 
data night. Since the HET is queue scheduled, stan- 
dards from subsequent nights were sometimes used. At 
all other telescopes, multiple standard stars were ob- 
served per night under varying atmospheric conditions 
and different air-masses. The sensitivity function was in- 
terpolated between standard star observations when the 
solution was found to vary significantly with time. 

For blue objects, broad-coverage spectrographs can 
suffer significant second order contamination. In par- 
ticular, the standard HET configuration using a Schott 
GG385 long-pass filter permitted second-order effects 
redward of 7700 A. The effect on object spectra were 
small, but for blue WD spectrophotometric standards, 
second order corrections were needed for accurate deter- 
mination of the sensitivity function. This c orrec tion term 
was constructed following iSzokolv et al.l (|2004f ) . In ad- 
dition, since BL Lac spectra are generally simple power 
laws, we used BL Lacs observed during these runs (Shaw 
et al, in prep.) to monitor second order contamination 
and residual errors in the sensitivity function. This re- 
sulted in excellent, stable response functions for the ma- 
jor data sets. 

Spectra were corrected for at mospheric extinction us - 
ing standard values. We followed rKrisciunas et al.l ()1987l ) 
for WMKO LRIS spectra, and used the mean KPNO ex- 
tinction table from IRAF for P200 DBSP spectra. Our 
NTT, VLT, and HET spectra do not extend into the UV 
and so suffer only minor atmospheric extinction. These 
spectra were also corrected using the KPNO extinction 
tables. We removed Galactic extinction us ing IRAF's de- 
redde ning function and the Schlegel maps (Schlegcl et al. 
Il998f l. We made no attempt to remove intrinsic redden- 
ing (ie. from the host galaxy). 

Telluric templates were generated from the standard 
star observations in each night, with separate templates 
for the oxygen and water line complexes. We cor- 
rected separately for the telluric absorptions of these two 
species. We found that most telluric features divided out 
well, with significant residuals only apparent in spectra 
with very high S/N. On the HET spectra, residual sec- 
ond order contamination prevented complete removal of 
the strong water band red-ward of 9000 A. 

When we had multiple epochs of these final cleaned, 
flux-calibrated spectra with the same instrumental con- 
figuration, we checked for strong continuum variation. 
Spectra with comparable fluxes were then combined into 
a single best spectrum, with individual epochs weighted 
by S/N. 

Due to variable slit losses and changing conditions be- 
tween object and standard star exposures, we estimated 
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Table 1 

Observing Configurations 



Telescope 


Instrument 


Resolution 


Slit Width 


Objects 


Filter 
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A 








A 


A 


HET 


LRS 


15 


2 


77 


GG385 


4150 


10500 


HET 


LRS 


8 


1 


1 


GG385 


4150 


10500 


McD 2.7m 


LCS 


15 


2 


1 




4200 


8200 


NTT 


EFOSC2 


16 


1 


14 




3400 


7400 


NTT 


EMMI 


12 


1 


8 




4000 


9300 


Palomar 200" 


DBSP 


5/15 


1 


3 




3100 


8100 


Palomar 200" 


DBSP 


5/15 


1.5 


1 




3100 


8100 


Palomar 200" 


DBSP 


5/9 


1.5 


2 




3100 


8100 


VLT 


FORS2 


11 


1 


22 




3400 


9600 


VLT 


FORS2 


17 


1.6 


2 




3400 


9600 


WMKO 


LRIS 


4/7 


1 


11 




3100 


10500 


WMKO 


LRIS 


4/9 


1 


22 




3100 


10500 



Note. — For DBSP and LRIS the blue and red channels are split by a dichroic at 
5600 A; the listed resolutions are for blue and red side, respectively. 



that the accuracy of our a bsolute spectrophotometry is 
~ 30% (jHealev et al.l [20081 ) . although the relative spec- 
trophotometry is considerably better. 

Redshifts were generally confirmed by multiple emis- 
sion lines and derived by cross-correlation analysis us- 
ing the rvsao package (jKurtz fc Minklll998D . For a few 
objects only single emission lines were measured with 
high S/N. In general we could use the lack of otherwise 
expected features to identify the species and the red- 
shift with high confidence. Nevertheless, single-line red- 
shifts are marked (Tables [3] and 0} with a colon (:) and 
discussed in ij2.4l Velocities are not corrected to helio- 
centric or LSR frames. 

Reduced spectra of the newly measured objects are 
presented in Figure [T] (full figure available on-line). 

Fig. Set 1. Spectra 

2.4. Individual Objects 

A few observations require individual comment: 
J0023+4456: This spectrum had low S/N (~ 2) and, 
while C IV, C III, and Mg II lines were tentatively iden- 
tified, we mark the final solution as uncertain. 

J0654+5042: This blazar, 0.5" from a bright star, was 
observed using DBSP with the slit along the blazar-star 
axis. We were able to isolate the stellar light by extrac- 
tion from the wings of the convolved object. This spec- 
trum was scaled and subtracted from the blazar spectrum 
to remove the stellar rest wavelength absorption features, 
giving a relatively clean blazar spectrum. While the 
emission line and redshift measurements are unambigu- 
ous, the continuum spectrophotometry should be treated 
with caution. 

J0949+1752, J1001+2911, J1043+2408, and 
J1058+0133: These blazars show only a single 
broad emission line; good spectral coverage gives 
high confidence of identification as Mg II. 

J1330+5202: This blazar shows only strong [O II] in 
emission at 3727 A, with associated host Ca H/K ab- 
sorptions. No broad emission lines are observed. 

J1357+7643: Only one broad line is detected with 
high significance. A weaker corroborating feature im- 
plies C IV, Mg II at z=1.585. However, a z=0.431 so- 
lution cannot be ruled out so we mark this solution as 
tentative. 

J2250-2806: With a single strong narrow line (assigned 



to [O II] 3727 A), this blazar has an uncertain redshift. 
No broad emission lines are observed. 

3. MEASURED PROPERTIES 

3.1. Primary Spectra and Multiple Epochs 

For some objects, multiple exposures at different 
epochs show significant spectral variability. In general, 
follow-up epochs have higher resolution or S/N, or tar- 
geted the blazar in a lower flux state. In a few cases SDSS 
spectra were published after we had obtained indepen- 
dent data. We adopt the spectrum showing the highest 
S/N for broad line detection as the primary spectrum 
used in this analysis. In general, this is the spectrum 
used to solve for the source redshift. As noted above, 
for particularly faint objects we can make S/N- weighted 
combinations of spectra (for which no strong flux varia- 
tion is seen) to obtain the best primary spectrum. This is 
particularly valuable for the HET data, as limited track 
time requires short observations, but queue scheduling 
with very stable configurations allows easy combination 
of multiple epochs. 

As noted above, in six cases, sources in this sample 
transitioned from a bright (continuum dominated) state 
to a lower state. For these the 'Primary' spectrum is from 
the low state, allowing better line measurements, during 
which the source is a nominal FSRQ (broad line EW >5 
A). We nevertheless retain the original BL Lac typing 
(see §3.7\i . Clearly a more physical separation between 
the classes, and attention to duty cycles in the various 
'states' is needed. 

3.2. Continuum Properties 

We have fit power law contin ua to all our spectr a, us- 
ing the scipy. optimize. leastsq (pones et al.l 12001— H rou- 
tine based on Levenberg-Marquardt fitting to estimate 
parameter values and statistical errors. For these fits 
we first excise regions around the strongest emission 
lines in the blazar rest frame: AA1200-1270, 1380-1420, 
1520-1580, 1830-1970, 2700-2900, 3700-3750, 4070- 
4130, 4300-4380, 4800-5050, 6500-6800 A. We then fit 
simple power law flux spectra to the entire spectral range, 
excluding regions at the blue or red end with unchar- 
acteristically large noise. We do not here attempt to 
separate the thermal disk contribution. The results are 
presented as power law indices, a, and continuum fluxes 
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Wavelength (A) 

Figure 1. Spectra of the Fermi blazars reported on in this work, 
All spectra available in electronic figure set. 

at 10 14,7 Hz (~ 5980 A), the center of our spectral range, 
as measured in the observer frame. These data may be 
combined with multi wavelength data to study the SED 
of the blazers in our sample. 

The statistical errors on F v (Tables [3] and are gen- 
erally small, but should be convolved with the overall 
fluxing errors, estimated at 30%. We estimate errors on 
the spectral index by independently fitting the red and 
blue halves of the measured spectrum. We then sum the 
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ordered by RA. Redshifts and optical types are listed for each spectrum. 

differences in quadrature for an estimated error: 

«err = (K- e d ~ a f + ("blue ~ a f) ^ W 

Note that large errors bars generally indicate deviation 
from power-law continua rather than poor statistics. We 
find that extreme values of a seem to correlate with large 
Galactic Ay , suggesting errors in the assumed extinction 
and residual curvature in the corrected spectrum. 

3.3. Line Properties 
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All line measurements are conducted in the object's 
rest frame. The local continuum and line parameters are 
measured with the scipy.optimize.leastsq routines. We 
first isolate the local continua surrounding the line, as in 
£13.21 (including an additional pseudo-continuum for Mg II 
from the Fe II/Fe III line complex). Then, the line is fit 
to the continuum-subtracted data. Reported errors are 
purely stati stical. 

We follow [Shen et all (|2011f) in their line measurement 
techniques to facilitate comparison of our results with 
with measurements in the extensive SDSS quasar cata- 
log. We focus on lines used for black hole mass estimates: 
H/3 at 4861 A, Mg II at 2800 A, and C IV at 1550 A. 

Before reporting kinematic widths, we subtract in 
quadrature the observational resolution, to eliminate bias 
between objects measured by different telescopes, and to 
improve estimates of narrower lines. In higher redshift 
sources, emission lines are sometimes contaminated by 
intervening or associated absorption line systems. We 
visually check all emission lines. 

We found strong absorption lines (associated and self- 
absorption systems) in a number of objects and edited 
these out before fitting the emission lines. Associated 

Mg II absorbers were found in J0043+3426 (Av 2400 

km/s), J0252-2219 (Av 3500 km/s and Av 2800 

km/s), J1120+0704 (Av 1600 km/s), J1639+4705 

(Av ~ -600 km/s) and J2212+2355 (Av ~ -2400 
km/s). In J0325+2224 we find self-absorbed C IV ab- 
sorption. In J2321+3204, we find two systems (Av ~ 
—400 km/s and Av ~ 100 km/s) and remove the absorp- 
tions from both Mg II and C IV before fitting. Finally, 
in J2139-6732, we edit out an intervening (z = 0.923) 
Fe II absorption system coincident with the C IV emis- 
sion line. We next summarize particular issues for the 
three major species fit. 

3.4. Fitting H(3 

For H/3, it is essential to include narrow com- 
ponents in the line fit. After power law contin- 
uum removal, we simultaneously fit broad and nar- 
ro w H/3, narrow [Q III]495 9 A, and narrow [O III] 5007 
A (|McLure fe Dunlopll2004h . We fix the rest wavelengths 
of narrow H/3 and the [O III] lines at the laboratory val- 
ues and fix widths at the spectral resolution, as measured 
from sky lines; the broad H/3 center and width are free to 
vary. All lines are modeled with single Gaussian profiles. 
The continuum is measured at 5100 A. 

3.5. Fitting Mg II 

For Mg II, we fit the pseudo- continuum from broad 
Fe II/Fe III alon g with a power law. For the former, we 
use a template (|Vestergaard fc Wilkesll200l convolved 
with the observational resolution. The power law is fit 
using the AA2200-2700 and 2900-3100 A portions of the 
spectrum; both continuum component fluxes and the 
power law index vary freely. 

After continuum subtraction, we fit the Mg II line 
with broad and narrow Gaussian components. For this 
line, the continuum measurement is made at 3000 A, to 
minim ize residual Fe contamination (|McLure fe Dunlod 
12001 . 

3.6. Fitting C IV 



Here, we fit the power law continuum over AA1445- 
1464 A and 1700-1705 A. After subtraction we fit the 
C IV line with three Gaussians, following iShen et aLl 
(|2011D . and report the full FWHM of the line. Any 
narrow self-absorption components are visually identi- 
fied and removed prior to the fit, as with the intervening 
absorbers described above. 

3.7. Objects with Special Type Classifications 

Several blazars were classified as BL Lacs in initial 
epoch observations. At the 'primary' spectrum epoch, 
with low continuum, each was a nominal FSRQ. The 
objects which changed (and continuum decrease) were: 
J0058+3311 (8x), J0923+4125 (4x), J1001+2911 (6x), 
J1607+1551 (5x), J2031+1219 (4x) and J2244+4057 
(lOx). 

With very high S/N observations, we were able to de- 
tect broad lines at high significance at EW levels < 5 A 
in several objects. These were thus 'BL Lacs' at all of 
our epochs, but can be analyzed along with the FSRQ. 
The BL Lacs (and strongest broad line EWs) were: 
J0430-2507 (Mg II at EW=0.9 A), J0516-6207 (C IV 
at EW=1.6 A; C III, Mg II also present), J1058+0133 
(Mg II at EW=2.2 A), J2236+2828 (Mg II at EW=4.9 A) 
and J2315-5018 (Mg II at EW=3.8 A). These EW mea- 
surements are in observed frame. Clearly these sources 
are transitional between our standard FSRQ and BL Lac 
types. 

3.8. Calibrating our Measurements 

For 53 of the 64 SDSS spectra in our sample lShen et all 
(|2011l) tabulate standard measurements of the continuum 
and emission lines properties. This allows us to check the 
consistency of our measurements techniques. 

We find that our measurements of conti nuum and 
line lu minosities are consistent with those of IShen et al.l 
(pOll . although some scattered differences in line lumi- 
nosity are seen, as the luminosity is sensitive to the fitting 
details. In particular, for Mg II and H/3, the variations in 
the narrow line component flux can change the measured 
broad line luminosity. Nevertheless, we see no system- 
atic o ffset between our results and those of IShen et al.l 

(pon . 

However, we do note that our fitted kinematic FWHM 
values are systematic ally ~ 10% large r (for the same 
sources) than those of IShen et al.l (|2011f) . For Mg II and 
H/3, we attribute this to small differences in the fit to the 
narrow component. For C IV, three Gaussians are used 
and we suspect that our central (narrowest) Gaussian is 
systematically weaker than that found in the SDSS fits. 
As we shall see, interesting differences from the SDSS 
quasars include narrower lines, on average, so we conser- 
vatively chose not to bias our fitter to force agreement. 
Nevertheless FWHM effects may be suspect at the 10% 
scale. 

3.9. Comparison to SDSS Non-Blazars 

To place our Fermi sample in the context of the broader 
QSO population, we compare our measurements to those 
of the l argely (> 90%) radio quiet SDSS DR7 quasar 
sample (|Shen et al.ll201ll) . We find significant offsets in 
the mean values of the continuum luminosities, line lu- 
minosities, and line FWHMs of the two classes as plotted 
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Figure 2. Line and continuum properties of our 1LAC sample (thick colored lines) compared to the SDSS clean sample (Shcn ct al. 2011) 
(thin black lines). To match the flux limit of the clean sample, we additionally plot the subset of Fermi objects with magnitude < 19.1 
(dotted lines). For each species and parameter, the plots are labeled with the SDSS mean value and the offset A and estimated statistical 
significance A/a of the means of the Fermi FSRQ, as measured in log space. The corresponding properties for the flux-limited sub-sample 
are given in parentheses. 
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in Figure [5J Since our sample extends fainter than the 
SDSS QSOs, we also plot the distributions for a magni- 
tude cut i < 19.1 sample (|Schneider et al.ll2010D . 

To quantify these differences, we compute the medians, 
as well as the semi interquartile range (SIQR) for each 
distribution. We can use these values to estimate sta- 
tistical errors on the median as a — SIQR/y^n, where 
n is the number of points in the sample. In Figure [2j 
we report both the median offset, A, and its estimated 
significance, A/er, for each parameter. The correspond- 
ing medians and errors for the magnitude-cut sample are 
given in parentheses. 

Several distributions show large offsets from the SDSS 
population. Our C IV and Mg II lines are significantly 
less luminous than those of the SDSS QSOs. Comparing 
the values for the flux-cut sample, we see that a signif- 
icant difference persists. We conclude that we are sam- 
pling a population with intrinsically weaker lines or that 
fainter, weaker-lines objects are lifted into our samples 
by the addition of extra non-thermal continuum flux. It 
seems that the local continua around Mg II and C IV are 
similar for our objects; we do notice that our H/3 con- 
tinua seem substantially brighter than those of the QSO 
population. We suspect that a lower redshift distribu- 
tion, and larger host-light contribution (see below) may 
explain this effect. 

For all species (and for the sub-sample) we find that 
our measured lines are significantly narrower than those 
the SDSS average. We attribute this to preferential align- 
ment, further discussed in S£SJ 

4. EXCESS CONTINUUM 

We expect our radio-7 selected population to have 
significant continuum contribution from non-thermal jet 
emission. At low redshift, there may also be additional 
galaxy light, especially around the H/3 line, to the red of 
the Balmer break. 

4.1. Estimated Continua from Lines 

In order to characterize the non-thermal emission in 
our spectra, we estimate a predicted continuum luminos- 
ity from the emission line fl uxes, scaling to a non-blazar 
sample. iShen et al.l (|2011l ) reports line and continuum 
luminosities for the DR7 quasars. These are highly cor- 
related (rxy = 0.85, 0.91, and 0.68 for H/3, Mg II, and 
C IV respectively), and they fit continuum luminosities 
to line luminosities for Mg II and C IV, finding 

log L 3000 = 1.016 ± 0.003 • log L Mg n + 1.22 ± 0.11 (2) 

log Li 350 = 0.863 ± 0.009 ■ log L c IV + 7.66 ± 0.41 (3) 

We extend this treatment by fitting the SDSS QSO re- 
sults for the H/3 line: 

log L 5100 = 0.802 ± 0.049 • log L Rfj + 1.574 ± 0.060 (4) 

Since the line flux continuum predictions are calibrated 
to a non-blazar sample, they should be independent of 
the non-thermal jet activity. Of course for both samples, 
at low redshift, there may still be galaxy light contam- 
inating the H/3 region, increasing the scatter. We com- 
pare the predicted to observed continuum levels for our 
Fermi FSRQ, finding that the lines predict (on average) 
44%, 40%, and 32% of observed continuum luminosities 
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Figure 3. The non-thermal dominance is the ratio of the mea- 
sured to the predicted continuum emission near each of the lines. 
The H/3 sample may be contaminated by galaxy luminosity as well. 
BL Lacs in the primary spectrum are circles. Squares indicate BL 
Lac designations from a previous spectrum. Large non-thermal 
dominance is weakly correlated to low spectral indices in Fermi 
observations. We thin error bars for NTD < 0.8 for readability. 

for H/3, Mg II, and C IV respectively. We attribute the 
excess to non-thermal 'jet' emission. 

It is of interest to inspect the non-thermal dominance 
(NTD = Lobs/ Lp re d) for individual objects in our sam- 
ple. While most objects have NTD ~ 1 (i.e. mostly 
thermal continua), for 57 blazars we have NTD > 2 and 
18 spectra show NTD > 10 based on measurements of 
at least one line. We find significant correlation in same- 
object NTD measured from the three species (rxy = 
0.51 for Mg II and C IV; rxy=0.73 for H/3 and Mg II) H/3 
generally measures larger NTD, which we attribute to 
contaminating galaxy light. Thus for 10% of our sample, 
the spectra are BL Lac-like and analysis of the optical 
properties must account for this dominant non-thermal 
emission. 

4.2. Non-thermal Continuum Pollution 

In Figure 02 we report on the NTD of our sample. An 
optica lly selected population, such as that of IShen et al.1 
(|201lL peaks at log NTD = 0, the line- predicted ther- 
mal flux is comparable to the measured flux. Unsurpris- 
ingly the Fermi blazars extend to much larger NTD lev- 
els. This plays an important role in continuum-calibrated 
estimates of hole mass in <|5] 

We further find a weak, but interesting correlation be- 
tween NTD and Fermi 7-ray spectral index. The Fermi 
FSRQ population has < 7 >= 2.44, whereas the BL Lac 
population has < 7 >= 2.04. Thus the 7-ray hardness is 
a good predictor of the continuum strength in the broad 
band SED. Indeed, we find that the objects in our sample 
with the largest NTD show relatively hard Fermi spec- 
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Figure 4. On the left, masses from the lines are plotted against traditional virial masses for Fermi FSRQ. Most objects lie to the right 
of the 1:1 (solid) line, indicating larger traditional (continuum) masses. The dashed line shows the effect of a 3x NTD. Closed symbols 
indicate objects with error < 0.4 dex in both masses. In the right panel, objects with multiple observations are plotted. For the three 
objects with significant (> 2a) mass offset, we connect the two masses with a solid line, and plot the error ellipse around the primary 
observation. For others, we measure the continuum luminosity in both spectra, and mark the mass offset from continuum change by a 
dotted line. 



tra. It will be of interest to trace this trend deeper into 
the BL Lac population, although NTD estimates will be 
more difficult for this sample. Note that while actual 
BL Lacs (circles), unsurprisingly, are mostly in the high 
NTD regime, a number of nominal FSRQ are found in 
this region, as well. 

5. THE CENTRAL ENGINE 

5.1. Traditional Virial Mass Estimates 

As in IShen et ahl (|2011f ) , we estimate traditional virial 
black hole masses from a relation of the form: 



log 



( M, 



BH 



blog(XL x ) + 21og(FWHM) (5) 



where a and b are calibrate d from reverberation map - 
ping for each line species (jMcLure fc Dunlod [2004). 
The broad line FWHM, in km s _1 , and the con- 
tinuum luminosity XL\, in units of 10 44 erg s _1 , 
are measured as described in ^3.41 i)3.51 and £)3.6I 
Values of a and b ar e tabu lated in Table [21 from 
IVestergaard fc Peterson! (1200(1 Therea fter VP06) for H 
3 and C IV and iVestergaard fc Osmerl (|2009fl (hereafter, 
VO09) for Mg II. A scatter of ~ 0.4 dex has b een inferred 
for vi rial masses in optically selected samples (IShen et al.l 

unni). 

For our Fermi sample we estimate mass from H 3 for 
50 objects, from Mg II for 176 objects, and from C IV 
for 68 objects. For 39 and 50 objects respectively, we 
measure mass from both H and Mg II, and from Mg II 
and C IV. Masses are reported in Table [3] and Table Si 

We urge caution in a naive interpretation of these 



masses, as applied to blazers. Due to the significant 
NTD of our sample, we expect continuum luminosity 
to not scale as in the reverberation mapping sample, as 
discussed below in §5.21 Further, preferential alignment 
and a non-spherical BLR yields systematically narrower 
lines, as discussed in SjSJ 

5.2. Mass Estimates from Lines Alone 

Bearing in mind the substantial NTD of much of 
our Fermi sample, we wish to make mass estimates 
from the line l umin osity an d line kinematic w idth alone. 
iGreene fc Hoi (|2005l ) and iKong et all (|200l calibrate 
such estim a tors f rom reverberation mapping samples. 
IShen et all (|2011[ ) calibrates estimators for Mg II and 
C IV from the SDSS sample. We have augmented this, 
using the SDSS mass estimators and developed our own 
coefficients for H/3 with values consistent with the pre- 
dicted continuum luminosity and mass estimates of ij5.ll 
The coefficients for these estimators are listed in Table [2] 
for a formula of the form of Equation [5J replacing XL\ 
with line luminosity in units of 10 44 erg s _1 . 

In Figure 2] we compare the line-estimated masses for 
our radio-7 selected sample with estimates from the tra- 
ditional continuum virial mass equations. We find the 
line masses are smaller than the traditional virial masses 
by 0.14 dex on average. This is a small, but significant 
effect. More interesting is the comparison for different 
species and different epochs. For example the nearby, low 
mass, low luminosity H/3 mass sample shows a larger av- 
erage line-mass decrease than the powerful C IV FSRQ. 
We attribute this primarily to host light pollution in the 
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Table 2 

Coefficients for Black Hole Mass Estimates 



Source 


H/3 


Mg II 


C IV 


Reference 




a b 


a b 


a b 




Mass from Continuum 


0.672 0.61 


0.505 0.62 


0.660 0.53 


MD04 (H/3,C IV), VO09 (Mg II) 


Mass from Lines 


1.63 ± .04 0.49 ± .03 


1.70 ±.07 .63 ±.00 


1.52 ±.22 .46 ±.01 


this work (KB). Shen et al. (20111 (Mg II. C IV) 



Note. — Coefficients for the vir ial fitting rel ation in Equation \5\ Coefficients derived from this work, IMcLure &: Dunlopl (|2004l . 



IMote . — Oo cmcicnts lor tne v irial nttmg relation m liquation | 
MD04). IVestergaard fc Osmeil flM VO09), and lShen et aJTlpOill ) 
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Figure 5. Blazar black hole (line) masses as a function of red- 
shift. Masses are over-plotted on the Sloan average quasar masses. 
Objects with i < 19.1 are plotted as open points. Large green 
crosses show the average Mgll-estimated mass and range for our 
magnitude-cut Fermi blazar sample in six redshift bins. 

former and thermal disk domination in the latter. 

For those objects where we have multiple epochs of ob- 
servation (see i j3.ll for details), we calculate the mass in 
both spectra. Due the statistical uncertainty in measur- 
ing the line kinematic width, and the large role that plays 
in mass calculations, in all but three cases, the measure- 
ment in the multiple spectra are not significantly differ- 
ent from that in the primary spectrum. In the right panel 
of Figure HI we also plot estimated mass changes from 
continuum luminosity fluctuations in spectra too poorly 
measured to estimate line properties. The general trend 
is for variations to extend above the line-determined val- 
ues, which show smaller fluctuation (i.e. continuum vari- 
ations dominate and bias mass estimates, as expected). 

5.3. Comparison to Optically Selected Quasars 

Using the less biased masses derived using line 
strengths, we compare in Figure [5] o ur radio-7 sampl e 
to the optically selected SDSS QSO (|Shen et all 120111 ). 
Interestingly, the mean mass varies with redshift in par- 
allel to the QSO sample, but lower. For the full Fermi 
sample the offset is 0.44 dex. For the magnitude-cut 
sub-sample (directly comparable to the SDSS QSO), the 
offset is 0.34 dex. 

The origin of this offset is unclear. Certainly the NTD 
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Figure 6. The Eddington ratio plotted against black hole mass. 
The green ellipse shows the centroid and approximate la range for 
values determined from Mg II lines for the SDSS QSO sample. The 
green line shows the approximate flux limit for our sample. 

of our sample causes some intrinsically fainter, less mas- 
sive blazars to be boosted into our flux cut sample. Fur- 
ther, our Fermi objects sample a preferred orientation 
(see . Another possibility is that these sources are on 
average in a more active accretion state than the typical 
QSO (see ^5.4[) . We also allow that the brighter blazars 
with with historical spectra in the literature may have 
systematically larger hole masses, decreasing the magni- 
tude of this offset. 

5.4. Eddington Ratio 

Following the unified model (|Urrv fc Padovarill [r995). 
we expect the quasi-isotropic thermal emissions of 
blazars to be similar to those in the underlying QSO 
population. 

To probe this, we measure the disc luminosity in Ed- 
dington units, LiEdd, using line measurements alone. We 
first use the line fluxes to estimate local continuum lumi- 
nosities, as above, t hen convert the s e to b olometric fluxes 
using factors from iRichards et all ((2006) (Lb i/L 5W Q = 
9.26, L bo i/L 3mo = 5.15, L bo i/Li 350 = 3.81). This pro- 
vides a bolometric thermal (= accretion power) luminos- 
ity, independent of any non-thermal contribution. This 
luminosity may be compared with the line-estimated hole 
masses to derive the Eddington ratio Fboil^Edd for the 
population. Figure [6] displays this ratio against line- 
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determined hole mass. Our sample has a systematically 
higher Eddington ratio and lower mass than the optically 
selected SDSS quasars. The radio-7 selection probes an 
intrinsically more active population than an optically se- 
le cted sample, both in f lux limited samples. 

iGhisellini et al.1 (|2011| ) has proposed that FSRQ are 
more active accretors than BL Lacs. While our high 
Eddington ratios for this FSRQ support this picture, we 
do not find any inverse correlation between Eddington 
ratio and non-thermal dominance as might be expected 
if high values indicate a transition to a BL Lac-type state. 

6. ORIENTATION OF THE BLR 

The virial black hole mass estimates require that the 
broad line region (BLR) has an isotropic velocity dis- 
tribution (S alviander et al.l l2007j ). In practice, there 
is a lack of consensus as to the shape of the BLR 
(jMcLure fc Dunlop! 12001 . With a more disc-like BLR, 
objects observed perpendicular to the disc will have 
lower kinematic FWHMs, and t hus, lower inferred BH 
masses. Following Dec arli et al.l (|2011f) . we can interpret 
the FWHM offset of §3.91 as a geometric effect and esti- 
mate the /-value, as a measure of the shape of the broad 
line region. 

If we assume that the underlying broad line shape dis- 
tribution is the same as that of the SDSS QSOs, and use 
our i < 19.1 flux-cut sample, we find that our Mg II and 
C IV lines are 0.14 dex narrower. Recalling that our line 
widths seem to be systematically slightly lower than the 
SDSS line estimates by 0.04 dex, the true offset may be 
as high as 0.18 dex. 

We assume we are probing a population with the same 
mass distribution as SDSS. Then, for a 0.14 dex offset our 
f- value is 1.38 /sDSS' wrnc h i s inconsistent with both a 
spherical BLR and a geometrically thin disk. In fact / 
can be related to the typical disk thickness ratio H/R: 



1 



< / >= 0.5 



— + sin~ 



-1/2, 



(G) 



where the average is taken over the observer inclination 
angle 9. Since there is good evidence that the Fermi FS- 
RQs are highly aligned (within 5°; Ajello et al, submitted 
to ApJ) by assuming that the SDSS QSO are viewed uni- 
formly (in sin 9) away from the obscuring equatorial torus 
(at 9 > 60°), we obtain an estimate for H/R ~ 0.4. Since 
the function is very flat at low 9, the result is insensitive 
to the precise degree of Fermi FSRQs alignment. In fact 
the required axis ratio only increases to 0.44 if the Fermi 
FSRQ are aligned within 2°. The decreased line width, 
and hence the decreased inferred hole mass are fully con- 
sistent with the preferred alignment of these blazars as 
expected in the Standard Model. Thus while we do not 
exclude underlying masses differences, our measurements 
are best explained by similar optically selected and radio- 
7 selected black hole mass distributions, but with viewing 
angles strongly aligned to the disk (and jet) axis for the 
latter. As noted in i)5.4l we do however find that the 
FSRQ show a higher Eddington ratio than the typical 
QSO, so some difference in accretion activity is indicated. 

7. COMPARISON WITH RADIO ACTIVITY 

Since the Fermi FSRQ associations are largely selected 
from a sample of radio bright, flat-spectrum sources, it 
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Figure 7. The black hole mass is uncorrelated to the intrinsic 

modulation in radio. 

may be useful to compare the optical spectroscopic prop- 
erties with radio measurements indicating strong non- 
thermal activity. One good measure of this activity is 
the radio variability. Much of the Dec > -20 blazar sam- 
ple has been monitored at 15 GHz at OVRO since before 
the launch of Fermi. These data have been used to mea- 
sure the intrinsic modulation index, a measure of the 
activity: 

m = -J- (7) 

Do 

where So and gq are the intrinsic mean flux density 
and its standard deviation of the radio light curve, 
estimate d from the data us ing a maximum-likelihood 
method ([Richards et al.ll2011[ ). It has, for example, been 
found that this index is smaller for high redshift FSRQ 
than lower redshift BL Lacs. 

We might, therefore, expect this modulation to cor- 
relate with either the black hole mass Figure [7] or the 
NTD Figure [8] The lowest intrinsic modulations do ap- 
pear among the highest mass holes, but there does not 
seem to be any increase in modulation for FSRQ with the 
higher NTD. Thus jet variability does not seem directly 
correlated with black hole and broad-line properties. It is 
possible that more significant differences may be traced 
to the radio variability timescale, but this requires longer 
radio monitoring and a more complete analysis of the 
flaring activity. 

8. CONCLUSIONS 

We have used the optical spectral properties of our 
measured Fermi FSRQ sample to characterize the state 
of the central engine. We find that the optical continuum 
is significantly augmented by non-thermal (synchrotron) 
emission, presumably associated with the jet. In fact, 
such emission dominates the optical continuum for a 
third of our measured FSRQ. Unsurprisingly the degree 
of NTD correlates with the 7-ray spectral index and the 
likelihood of BL Lac classification. 

To avoid bias from the additional continuum, we have 
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Figure 8. The intrinsic modulation is uncorrelated to the non- 
thermal dominance of the optical spectrum. 

developed estimates of the black hole mass and the Ed- 
dington luminosity ratio which only use the line proper- 
ties. These measurements show that the Fermi FSRQ are 
more active accretors than the bulk QSO population (at 
a given mass). We also find evidence for a significantly 
lower mean measured BH mass, even for a sample with a 
similar flux limit. Intriguingly, we find that the redshift 
evolution of the mass tracks that of the full QSO pop- 
ulation. Thus, although non-thermal flux can pull less 
luminous, lower mass AGN into the sample, we specu- 
late that the primary effect is due to the high expected 
degree of alignment for the FSRQ sample. If alignment 
is indeed the explanation, the data suggest a modest de- 
gree (H/R ~ 0.4) of flattening in the broad line region, 
and that we probe a similar mean mass population (and 
evolution) as the bulk QSO distribution. 

While the present sample shows no strong correla- 
tion with the bulk radio variability, we expect that a 
broader look at the Fermi blazar sample, including com- 
parison with the BL Lac properties and study of the BL 
Lac state-duty cycle and timescales for radio and optical 
modulation, will reveal additional correlations with the 
properties of the central black hole and its surrounding 
broad line region. This extension to the BL Lac portion 
of the sample is under way. 
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of the University of Texas at Austin, the Pennsyl- 
vania State University, Stanford University, Ludwig- 
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E. Eberly. The Marcario Low Resolution Spectrograph is 
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Table 3 

Line Properties for z < 1 Broad Line 1LAC AGN 



H/3 



Mg II 





io- 


28 erg cm _2 s 


- 1 Hz _1 


-^5100 

erg s -1 


erg s _1 


kms — 1 




-^3000 

erg s — 1 


-^Mg II 
erg s 


i_ VV 11 IV J. 

kms 


M ass 

IV J. cloo 


_L <^ 1 o vJ 1} \3 


J0004-4736 


0.880 


11.4±0.0 


-0.99±0.08 










45.339±0.011 


42.896±0.151 


2700±500 


7.85±0.36 


NTT 


J0008+1450 


0.045 


73.8±0.1 


-1.79±0.23 


43.216±0.002 


40.921±0.120 


10800±2600 


8.19±0.48 










SDSS 


J0017-0512 


0.226 


30.4±0.1 


-0.46±0.36 


44.353±0.004 


42.388±0.107 


2300±500 


7.55±0.45 










HET 


J0024+0349 


0.545 


27.1±0.0 


-1.88±0.08 










45.110±0.001 


42.585±0.057 


3000±200 


7.76±0.14 


WMKO 


J0043+3426 


0.966 


6.6±0.0 


-1.24±0.03 










45.201±0.006 


42.807±0.072 


3400±300 


8.01±0.16 


WMKO 


J0050-0452 


0.922 


7.7±0.1 


-0.92±0.19 










45.198±0.013 


43.140±0.141 


3300±600 


8.20±0.34 


HET 


J0102+4214 


0.874 


30.1±0.1 


-1.63±0.86 


46.001±0.004 


43.460±0.038 


1900±200 


7.92±0.16 


45.409±0.007 


43.593±0.072 


3300±300 


8.49±0.17 


WMKO 


J0102+5824 


0.644 


167.0±0.3 


-1.25±0.36 










46.093±0.076 


43.449±0.256 


4100±1200 


8.57±0.61 


HET 


J0112+2244 


0.265 


18.6±0.2 


-2.71±0.79 


















HET 


J0113+1324 


0.685 


15.2±0.0 


-0.01±0.26 


45.063±0.006 


43.094±0.062 


3800±500 


8.35±0.25 


45.263±0.006 


43.389±0.082 


4300±400 


8.57±0.19 


SDSS 


Note. — Tabic 


[3]is publ 


ishcd in its cntii 


■ety in the electro 


nic edition of this . 


journal; A portion 


is shown here for gu 


idancc regarding its form and conti 











Table 4 

Line Properties for z > 1 Broad Line 1LAC AGN 



Name 


z 

10~ 


-28 er g cm — 2g 


a 

-iHz" 1 




Mg II 








C IV 






Telescope 


-^5100 

erg s 


L Mg II 

erg s 


EWHM 

kms -1 


Mass 


-^3000 

erg s 


L C IV 
erg s — 1 


EWHM 

kms -1 


Mass 


J0011+0057 


1.493 


5.2±0.1 


-0.51±0.05 


45.486±0.021 


43.541±0.081 


4900±500 


8.80±0.19 


45.800±0.069 


43.493±0.435 


2500±1300 


8.09±1.02 


SDSS 


J0015+1700 


1.709 


8.0±0.0 


-2.75±0.38 


46.070±0.004 


44.178±0.043 


5900±300 


9.36±0.10 


45.818±0.016 


44.180±0.068 


5900±500 


9.15±0.17 


WMKO 


J0023+4456 


2.023 : 


1.3±0.1 


-0.20±0.50 










45.213±1.000 


43.336±0.875 


1900±2200 


7.78±2.29 


HET 


J0042+2320 


1.426 


5.1±0.1 


-1.40±0.08 


45.533±0.017 


43.412±0.137 


6900±1100 


9.01±0.32 










HET 


J0044-8422 


1.032 


9.7±0.0 


-0.14±0.03 


45.401±0.007 


43.667±0.113 


3900±500 


8.68±0.26 










VLT 


J0048+2235 


1.161 


11.5±0.0 


-1.34±0.18 


45.639±0.002 


43.143±0.049 


4300±200 


8.43±0.11 


45.916±0.027 


43.272±0.115 


3400±500 


8.25±0.28 


WMKO 


J0058+3311 


1.369 


3.1±0.1 


-2.49±0.37 


45.331±0.003 


42.805±0.054 


3400±200 


8.01±0.13 


45.177±0.035 


43.451±0.076 


2200±200 


7.97±0.19 


WMKO 


J0104-2416 


1.747 


8.8±0.0 


-0.81±0.22 


45.817±0.010 


43.613±0.111 


5000±600 


8.85±0.26 


45.989±0.025 


44.121±0.292 


4900±1700 


8.97±0.69 


NTT 


J0157-4614 


2.287 


1.9±0.0 


-1.07±0.33 


45.632±0.030 


43.216±0.171 


2500±500 


7.98±0.40 


45.522±1.000 


44.083±7.015 


3000±24900 


8.52±20.15 


VLT 


J0226+0937 


2.605 


16.9±0.0 


-0.58±0.26 










46.811±0.006 


44.582±0.720 


8500±7500 


9.65±1.76 


P200 



CO 
CO 

% 



Note. — Table [4l is published in its entirety in the electronic edition of this journal; A portion is shown here for guidance regarding its form and content. 



